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Abstract 

We calculate open charm production in Pb+Pb reactions at SPS energies within 
the HSD transport approach - which is based on string, quark, diquark (q, q, qq, qq) 
and hadronic degrees of freedom - including the production of open charm pairs from 
secondary 'meson'- 'baryon' (or quark-diquark and antiquark-diquark) collisions. It 
is argued that at collision energies close to the cc pair threshold the dominant 
production mechanism is related to the two body (or quasi two body) reactions 
ttN — > D(D*)A C , (S c ). Estimates within the framework of the Quark-Gluon String 
model suggest cross sections of a few \xb for irN — > DA C in the region of 1 GeV above 
threshold. The dynamical transport calculations for Pb + Pb at 160 A- GeV indicate 
that the open charm enhancement reported by the NA50 Collaboration might be 
due to such secondary reaction mechanisms. 
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1 Introduction 



In the last decade the interest in hadronic states with charm flavors (c, c) has been rising 
continuously in line with the development of new experimental facilities M. This relates 
to the charm production cross section in pN and ttN reactions as well as to their inter- 
actions with baryons and mesons which determine their properties (spectral functions) 
in the hadronic medium. The charm quark degrees of freedom play an important role 
especially in the context of a phase transition to the quark-gluon plasma (QGP) j2| where 
charmonium (cc) states should no longer be formed due to color screening |3|, M. However, 
the suppression of J/ty and ty' mesons in the high density phase of nucleus-nucleus colli- 
sions U might also be attributed to inelastic comover scattering (cf. J7], ||, ||, [U], [LJJ and 
Refs. therein) provided that the corresponding J/\l/-hadron cross sections are in the order 
of a few mb [|1^, 13, |14], 15, HJ. Present theoretical estimates here differ by more than 
an order of magnitude j[7| especially with respect to J/^-meson scattering such that the 
question of charmonium suppression is not yet settled. On the other hand, the enhance- 
ment of 'intermediate-mass dileptons' in Pb+Pb collisions at the SPS has been tentatively 
attributed to an enhancement of 'open charm' in nucleus-nucleus collisions relative to pA 
reactions at the same invariant energy yfs [18 . Thus 'charmonium suppression' and 'open 



charm enhancement' are present facets of relativistic heavy-ion collisions, which provide 
a theoretical |L5|, |2(], [21], |22], |23| and experimental challenge for the future J23], [25"| . 

In this letter we will argue that the open charm excess reported in Ref. |18| might 
be due to secondary meson-baryon (quark-diquark) interactions. We will employ the 
HSD transport approach JF], [10|, [2j| for the overall reaction dynamics of nucleus-nucleus 



collisions using parametrizations for the elementary production channels for the charmed 
hadrons D,D,D*,D*,D S ,D S ,D* S ,D* S , J/^,^(2S),X2c from NN and irN collisions. We 
recall that in the HSD transport approach the initial stages of a irN, pp, pA or AA reaction 
(at high energy) are described by the excitation of color neutral strings, where the leading 
quarks and 'diquarks' in a baryonic string (or quarks and antiquarks in a mesonic string 
etc.) are allowed to rescatter again (in case of nuclear targets) with hadronic cross sections 
divided by the number of constituent quarks and antiquarks in the hadrons, respectively 



In order to explore the kinematical situation for secondary 'meson'-'baryon' interac- 
tions we show in Fig. 1 the distribution in the invariant collision energy \fs from such 
secondary interactions (solid histogram) for a central collision of Pb + Pb at 160 A-GeV 
from the HSD approach. This distribution extends to about 15 GeV while the correspond- 
ing primary baryon-baryon collisions are centered around 17.3 GeV (dashed histogram) 
and become less frequent than the 'meson'-'baryon' collisions for y/s < 13 GeV. Since the 
threshold for open charm production in meson-baryon collisions is ^/sq = + m^ c « 
4.15 GeV there are a lot of secondary collisions above threshold. The crucial question 
now is the magnitude of the open charm cross section in 'meson'-'baryon interactions for 
4.15 GeV < y/s < 15 GeV, i.e from threshold up to roughly 10 GeV above threshold. 

As argued in Refs. |20| , 2"8]| the creation of a cc pair at high \/s is due to a hard process 
and dominated by gluon-gluon fusion (cf. Fig. 2a) for an illustration of a pp reaction). On 
the other hand, the quark annihilation mechanism is found to contribute significantly or 



even to become dominant for irN reactions at lower \/s [29]. This mechanism is depicted 
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Table 1: The parameters a x , a and j3 in Eq. (p]) for inclusive D- meson production in nN 
reactions 



nN 


Meson 


[GeV] 


a x [mb] 


a 


P 


D» 


4.667 


0.273 


2.86 


1.28 


D° 


4.150 


0.247 


3.80 


1.26 


D+ 


4.671 


0.255 


3.22 


1.28 


D~ 


4.154 


0.286 


3.50 


1.22 


D°* 


4.951 


1.076 


3.14 


1.22 


D°* 


4.292 


0.774 


3.80 


1.26 


D + * 


4.955 


0.719 


2.86 


1.32 


D~* 


4.296 


0.839 


3.40 


1.24 


Df 


4.875 


0.0932 


3.62 


1.22 


D- 


4.435 


0.0545 


3.70 


1.34 


D+* 


5.162 


0.284 


3.42 


1.24 


D-* 


4.578 


0.163 


3.64 


1.34 



in Fig. 2b) for the particular final state DA C , where - by convention - D stands for a D- 
meson with a c-quark. However, the perturbative qq — > 1 gluon — > cc process is included 
in standard PYTHIA calculations |31J and drops very fast with decreasing y/s close to 
the threshold of D— meson production p8| . We recall that in Ref. [28| the open charm 



cross sections for pN and ttN reactions have been calculated within PYTHIA using MRS 



G (next to leading order) structure functions from the PDFLIB package |H] for the gluon 
distribution of the proton, a bare charm quark mass m c =1.5 GeV and — 1 GeV. The 
results from PYTHIA (scaled to the available data M) for D, D, D* D* D„ L>„ D* D* 



as a function of ^fs > 10 GeV have been parametrized by an expression of the form, 

cr x (s) = a x (l - Z) a Z-P, (1) 



s where 



denotes the threshold for the channel X in pN or ttN 



The individual parameters ax , ot and (3 from this fit are given in Table 1 

|28| the contribution to open charm production 
employing the parametrizations ([!]) - are about ~ 9 % in 



with Z = 
reactions. 

for 7iN reactions. As found in Ref. 
from 'secondary' interactions 
central Au + Au collisions at 160 A- GeV. This order of magnitude for the open charm 
enchancement is far below the observation of Ref. \TE\. 

On the other hand, close to threshold, one would expect the exclusive reaction ttN — > 
DA C to dominate. In case of 5* -wave production this cross section then should rise as 
~ (1 — Z) 0,5 for a two-body final channel. Such a behaviour is well known experimentally 
1321 for the reactions ir + p — > p + p, n~p — > un, -n~p — > K°A or n~p — > <pn as displayed in 
Fig. 3. These experimental cross sections can be described by the expression 



^7tN~>xb{Vs) = 



0.5/ 



7x 



(2) 
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with a p = 15 mb, a u = 12 mb, a^o = 2.5 mb, = 0.2 mb and 7^ = 7^ = 6, 7 P = 5 and 
70 = 10, respectively. In Eq. (||) y/so denotes the threshold for each reaction separately. 
For y/s — t/so > m n three-body final states (with an additional pion) become possible 
while multi-particle production dominates at high invariant energy. To demonstrate the 
relative importance of multi-particle final states the inclusive cross sections for p + ,u> and 
<p are shown in terms of the thick lines in Fig. 3 using the parametrizations from the 
review J7[. The inclusive yield for all open charm mesons with a c quark is shown in 
terms of the lower thick solid line in Fig. 3 where the parametrisations ([]]) have been 
employed with the parameters from Table 1 (cf. PSfl). The related data (full squares) 
have been taken from the review of Tavernier |33|. We recall that for — > 6 GeV 
the parametrization reflects the PYTHIA results (scaled to the available data) and at 
lower energies lacks a more fundamental justification as pointed out in Ref. |28|] . 

Since the binary exclusive reactions dominate the threshold behaviour for p, u, K° and 
4> production, one can expect a similar relation to hold also for the channels DA C , D*A C , 
Z)S C , D*£ c , etc. Here the reaction mechanism is expected to be dominated by D*-meson 
or .D*-Reggeon exchange as illustrated in Fig. 2c). The question that remains is the size 
of a related parameter ao in (Q). Here we refer to the Regge model for an estimate. 

In Ref. P|fl the cross section for the reaction ttN — > D(D*)A C has been estimated 



within the framework of the Quark-Gluon String Model (QGSM) developed in Ref. |35| . 
The GGSM is a nonperturbative approach based on a topological 1/Nf expansion in QCD, 
where Nf denotes the number of flavors, and on the color-tube model. This approach can 
be considered as a microscopic model describing Regge phenomenology in terms of quark 
degrees of freedom. It provides the possibility to establish relations between many soft 
hadronic reactions as well as masses and partial widths of resonances with different quark 
contents (see e.g. the review ]3Tj|). Recently, this model also has been successfully applied 



to the description of the nucleon and pion electromagnetic form factors |37]] and deuteron 
photodisintegration 



The amplitude of the reaction icN — > D(D*)A C corresponding to the planar graph of 
Fig. 2c) with u and c quark exchange in the t channel - using the Mandelstam variables 
s and t - can be written as (see Ref. |3l 



4r-|~J>A«(M) = V^glF^is/sfT^-^s/s), (3) 

where g ~ 5.8 is a universal flavor independent coupling constant, a u£ (t) = a is the 
D* Regge trajectory, s — 1 GeV 2 , Sq c = 4.9 GeV 2 is the flavor dependent scale factor and 
F(t) is the form factor describing the t dependence of the residue. We have considered 
two forms of the D* Regge trajectory, i.e. 

i) a linear trajectory (set a from Ref. 

a D *(t) = -0.86 + 0.5t (4) 

ii) a nonlinear trajectory, i.e. the square root trajectory from Ref. |39fl 

a D * (t) = a D * (0) - 7 [VT - VT - t] (5) 

with ao*(0) = —1.02, 7 = 3.65 GeV -1 , VT =3.91 GeV. Both choices lead to very similar 
total cross sections such that we do not discuss the two choices separately. The form 
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factor F(t) determining the t-dependence of the residue has been parametrized in Ref. 
33 as 



F(t)=T(l-a D *(t)) (6) 

and is motivated by dual models. The form factor (|6|) is convenient for an analytical 
continuation of the amplitude to the region of positive t, where the T function in (^) is 
an exponentially decreasing function of t. In particular, the form (||) was successfully 
used to relate the normalization of different planar graphs at t= with the widths of the 



mesons lying on the corresponding Regge trajectories. For example, in Ref. |4(J a width 
T(D* + -»■ D°n + + D + 7i°) = (3/2)T(D* + -»■ D°n + )= 20 keV was found. This value is not 



very different from the result based on the QCD sum rules |4T| giving ~ 48 keV, while 
the experimental upper limit is 131 keV. 

However, in the region of negative t the parametrization (|6]) exhibits a factorial growth 
(which is faster than exponential) and therefore is not acceptable. For t < we will use 
the convential parametrization 

F(t) = T(l - a D *(0)) exp(i? 2 t) (7) 

with R 2 = -T- 0.2 GeV~ 2 . 

The differential cross section for the reaction ir^p — ► DA C then is 

da n - p ^ DAc _ 1 1 _ 2 

where p c ™ denotes the pion momentum in the cms, and the total cross section is obtained 
by integrating (|8]) over the kinematical allowed regime. 

The result for the channel ir~~p — ► DA C (using R 2 = 0) is shown in Fig. 3 by the lower 
solid line as a function of yfs — -y/so, which can be well described by the function (0) using 
an = 0.027 mb and jd = 7w = Ik = 6 (lower dashed line). The ratio of the uj cross 
section to the 'estimated' D cross section close to threshold ( ^Js — ^JJq ~ 0.4 GeV) is ~ 
500; a similar ratio holds for the inclusive cross sections at \fs — y^s^ ~ 30 GeV where 
explicit data are available for the ttN — > DD + X reaction. Thus our estimate in the 
Regge model sounds reasonable, however, has to be controlled by experiment in order to 
allow for definite conclusions. 

We mention that similar exclusive cross sections can also be obtained from boson- 
exchange models with D* or D\ exchange (cf. Fig. 2c) employing monopole form factors 
at the vertices with cutoffs A (w 1.6 — 1.9 GeV) while fixing the couplings g^D* an d 
gD* 7r= D 1 to the Ti D and nD* decay widths, respectively. Furthermore, the cross section 
for the reaction iViV — > D(D*)A C N has been calculated within the Reggezied one-pion 
exchange model using the pion form factor F n = exp(R 2 t) with R 2 = 0.l(GeV/c)~ 2 as 
in Ref. [J53J]. Since the explicit cross sections for the iViV reactions at low ^/s are of no 
major importance for our present study, we refer the reader to a forthcoming publication 



on these issues 42 



For the transport calculations (to be discussed below) we estimate the cross section for 
charmed hadron production in ttN reactions by employing the ansatz @ for all charmed 
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hadrons such as D, D*, D s , D* with associated charmed hyperons A c , S c , A^.... We explic- 
itly adopt a D = l/3a D * = cld* = 3qd s (as in Ref. f28|) with Od ~ 0.027 mb. These 



estimates, of cause, only have exploratory character as stated above. 

We now turn to the results of transport calculations for the reaction Pb + Pb at 
160 A-GeV. A description of the transport approach is given in Refs. [|7j, |27| and the 
production and propagation of charmed mesons is described in Refs. |28], [|3[ The 
novel phenomenon addressed here with respect to Ref. is the secondary production 
of open charm mesons by 'meson'-'baryon' collisions with the new parametrizations (|2]) 
for the low-energy ttN — > Z)(Z)*)A C (S C ) processes. We mention in passing that the role of 
secondary reactions on intermediate mass dilepton pairs via the qq — > 7 — > l + l~ (Drell- 
Yan) mechanism has been investigated in the UrQMD transport model in Ref. ||5| for 
nucleus-nucleus collisions at SPS energies before. 

In Fig. 4 we show the open charm multiplicity (all D, D,D*, D*,D* S , D*,D S ,D S ) 
as a function of the impact parameter b from the HSD approach for Pb + Pb at 160 
A-GeV. Here the dashed line gives the yield from primary baryon-baryon collisions, which 
is essentially the yield from pn reactions at ^fs ~ 17 GeV times the number of 'hard' 



pN collisions (as given by Eq. (6) of Ref. [28]), while the solid line corresponds to the 



yield of open charm from secondary reactions as described above. Whereas for peripheral 
collisions the secondary production is slightly smaller than the direct channel, the yield 
from secondary meson-baryon channels becomes larger with decreasing impact parameter 
b. The tiny kinks in the curves in Fig. 4 are due to the limited statistics of the transport 
calculations. Without explicit representation we note that the final differential spectra 
of the .D-mesons in transverse momentum and rapidity are very similar to those from 
primary and secondary channels, respectively, due to rescattering with the surrounding 
hadrons (cf. Ref. 



In order to compare to the data from the NA50 Collaboration ][y| we have to adopt 
a model to convert the number of participating nucleons A part to the impact parameter b 
from the transport calculation. We use 

A part = 2A-N (b), (9) 

where N (b) stands for the number of noninteracting nucleons (with no 'hard' collisions) 
from the transport code at impact parameter b while A is the target (projectile) mass 
number. Denoting the expected number of D, D mesons from primary BB collisions by 
Nprim. and the number from secondary 'meson'-'baryon' interactions by N sec . this leads 
to the ratio 

R(b) = ^ + ^ (b) = R(A part (b)) = R( Apart), (10) 

which is displayed in Fig. 5 (solid histogram) in comparison to the data on the open 
charm enhancement from Ref. fl8|| . The calculated ratio first increases fast with A part 
and becomes almost constant for A part > 100. This is due to the fact that secondary 
semihard 'meson'-'baryon' reactions with 'wounded' nucleons, that have scattered at least 
once, already set in for large impact parameter b; their relative number then increases 
only slightly faster than ~ A part . Obviously, the general trend of the data can be roughly 
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described with increasing centrality for the secondary cross sections specified above. The 
present statistical uncertainty of the data [18|] does not allow for a final conclusion. 

In summary: in this letter we have argued that for 'low' energy 'meson'-'baryon' 
reactions the dominant cc production is related to the two-body (or quasi two-body) 
reactions ttN — > D(D*) A C (S C ). Estimates within the framework of the Quark-Gluon 
String model suggest cross sections of a few fib for ttN — » DA C in the region of about 1 
GeV above threshold. The estimated order of magnitude for the open charm cross section 
(cf. lower part of Fig. 3) is found to be compatible with the 'open charm enhancement' 
claimed by the NA50 Collaboration at the SPS |18[ without employing the assumption 



of thermal and chemical or statistical equilibrium as advocated in Refs. [|6|, |47], It 
should be stressed, however, that experimental investigations on open charm production 
in ttN reactions at invariant energies of 4.2 < s/s < 15 GeV are mandatory to confirm or 
disprove our suggestion. 

We, furthermore, note that the cross section for charmonia such as J/^f, Xc or ty' are 
not substantially enhanced by such secondary reaction channels since their cross sections 
are small in ttN collisions || such that no substantial 'enhancement' of charmonia 
relative to the primary NN reaction channels is expected. 



The authors are grateful to E.L. Bratkovskaya, A.B. Kaidalov and R. Vogt for clarifying 
discussions. 
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Figure 1: The distribution in the invariant collision energy y/s from primary baryon- 
baryon (dashed histogram; BB) and secondary 'meson'- 'baryon' interactions (solid his- 
togram; mB) for a central collision of Ph + Pb at 160 A-GeV from the HSD transport 
approach. The arrow at ^/s ~ 4.15 GeV denotes the threshold for D + A c production in 
secondary interactions. Note, that contrary to Fig. 13 in Ref. [|Hj 'meson' interactions 



with 'wounded' nucleons (or diquarks) are taken into account, too. 

10 



71 



N 



71 



N 



a) 




b) 



D 



c) 



_ < 



D 



A 



D 



A 



Figure 2: Diagrams for cc production in pp collisions via two-gluon fusion (a), for DA C 
production in ttN collisions via qq annihilation (b) and for DA C production in irN collisions 
via .D*-exchange (c). 
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Figure 3: The exclusive cross sections for the reactions n + p 

~p -> K° + A and n 
(0) (dashed lines) 



p + + p, n p — > + n. 



32| in comparison to the parametrizations 
The inclusive cross sections for p + 



it i> — j\ " ;\ and ~ +p — > + from Ref. 

, uj and cf) mesons are displayed in 
terms of the upper thick solid lines, respectively, within the parametrizations from the 
review [^]. The lower solid line is the prediction for the process 7r _ +p — ► D° + A c within 
the Regge model; the lower dashed line is the corresponding parametrization by (0). The 
inclusive yield for all open charm mesons with a c quark (D incl.) is shown in terms of 
the lower thick solid line where the functions (|I|) have been employed with the parameters 
from Table 1. The related inclusive data (*) have been taken from Ref. [33]. 
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b [fm] 



Figure 4: The multiplicity of all D, D mesons from the HSD transport approach for Pb+Pb 
at 160 A-GeV as a function of impact parameter b. The dashed line stands for the yield 
from primary baryon-baryon collisions while the solid line denotes the contribution from 
secondary 'meson'-'baryon' collisions within the parametrization (0). 
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Figure 5: The enhancement factor (|lCf ) for the production of all D, D mesons from the 
HSD transport approach for Pb + Pb at 160 A-GeV as a function of the number of 
participants A part (see text) in comparison to the data from Ref. [18 
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